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Abstract 

Recent progress in understanding the strong physics regime of QCD is described. 
The role played by condensates, particularly qq, in breaking chiral symmetry and 
generating constituent masses for u and d quarks is reviewed. The influence this has 
on hadrons with vacuum quantum numbers is emphasised. What we know of this 
sector from recent data on i;^)-radiative decays and from D decays to light hadrons 
is discussed. The key to further understanding is comprehensive analyses of such 
data, including that planned for Hall D at Jefferson Laboratory. 
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1 QCD vacuum 



To learn about the underlying theory of quark and gluon interactions we have to study 
hadronic collisions. In the debris of these lies the beauty of QCD. The task required to 
reveal this world of quarks and gluons is akin to an archaeological dig through the debris 
of Iraq in 2003 to unearth the civilisation of Babylon beneath. Over 30 years of digging 
considerable progress has been made in uncovering the strong physics aspects of QCD. 

What makes QCD more interesting than QED is the nature of the vacuum. For QED 
the vacuum is described by perturbation theory. It is essentially empty with a low density 
of particle-antiparticle pairs. In QCD because the interactions are stronger not only is the 
vacuum a denser sea of qq pairs and cloud of gluons, but so strong are the forces that 
condensates of quarks, antiquarks and gluons form in all colour singlet combinations. That 
of lowest dimension is the qq condensate. The scale of this condensate characterises the key 
non-perturbative effects of light hadron systems. Being non-zero it dynamically breaks the 
chiral symmetry of QCD. This simultaneously ensures that the pion is the Goldstone boson 
of this symmetry breaking and that the corresponding scalar field is the Higgs sector of the 
strong interaction, responsible for the masses of all light hadrons. 

The value of this condensate {uu) ^ {dd) can be found in at least three different 
ways m 12] , which here will be described as "phenomenological" , "experimental" and "the- 
oretical". Remarkably the 3 distinct ways give consistent results. The first and oldest 
phenomenological method is the application of the QCD sum-rules of Shifman, Vainshtein 
and Zakharov jHj to scalar and pseudoscalar currents. The sum-rules relate the matrix ele- 
ment of current correlators evaluated at low energies from hadronic data to their calculation 
at higher energies using the Operator Product Expansion. It is in this expansion that con- 
densates arise. Though such sum-rules have been studied for 25 years, recent precision has 
come from a better understanding of the use of contour improvement, of pinched weights in 
finite energy sum-rules and technological advances in the calculation of higher order correc- 
tions in perturbative QCD (see citations in Ref. 2). Agreement between the theoretical and 
experimental sides of the sum-rules gives {qq) ~ —(250 ± 25 MeV)^ at a scale of 2 GeV. Of 
course, such sum- rule analyses are only a consistency check on the size of condensates not 
an absolute determination. 

Much more direct experimental confirmation is obtained by measuring low energy tttt 
scattering precisely [21 El El E]- Pions being the Goldstone bosons of chiral symmetry breaking 
know about the size of {qq). This is reflected in the value of the vrvr scattering amplitude 
at the symmetry point in the middle of the Mandelstam triangle s = t = u = 4m^/3. At 
this unphysical kinematic point, the size of the amplitude increases by up to a factor of 4 as 
the qq condensate decreases in scale from 250 MeV ;5]- Precision measurements (to better 
than 10%) of vrvr interactions below cm. energies of 450 MeV have now become possible. 
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Combining such data with dispersion relations that incorporate the important constraint of 
the 3-channel crossing symmetry of tttt scattering allows the value of the amplitude at the 
symmetry point to be determined and hence fixes the size of the condensate While we 
await the precision measurement of the amplitude at threshold deduced from the lifetime 
of pionium [7j, one can use the difference of the phase of S and P-wave vrvr interactions as 
measured in K ^ ez/vrvr decay, as described in detail in Ref. 4. Results from the BNL-E852 
experiment [H] yield a condensate of ~ —(270 MeV)'^ at 2 GeV, showing that more than 90% 
of the Gell-Mann-Oakes-Renner relation for expanded in powers of the current quark 

mass is given by just the first term linear in rriq [H]. 

For a theoretical determination one can solve strong QCD in the continuum using the 
Schwinger-Dyson equations under certain plausible assumptions, as discussed in Refs. 10,1,2. 
These studies enable the behaviour of gluon, ghost and quark propagators and their interac- 
tions to be investigated in the small quark mass (or chiral) limit. Here considerable progress 
has been made in the past decade (TTHTU]. We understand that if the effective quark-gluon 
coupling becomes of order unity for momenta below 500 MeV or so, chiral symmetry is 
broken. Then a massless current quark (or one of mass of a few MeV), that propagates 
almost freely over very short distances, has an effective mass of 350 MeV at distance scales 
of 1 fm. The strong QCD dressing of the u and d quark propagators turns a current quark 
into a constituent quark. Remarkably, this behaviour corresponds in the chiral limit to 
the (qq) condensate having a scale of 250-270 MeV [12] , in reassuring consistency with the 
phenomenology just discussed. 

As emphasised repeatedly by Roberts and collaborators P^, the axial Ward identity 
ensures that the qq bound state with pseudoscalar quantum numbers is the Goldstone boson 
with its interactions governed by PCAC. In contrast the bound states with scalar and vec- 
tor quantum numbers have masses reflecting the mass of the fully dressed (or constituent) 
quark. The behaviour of the gluon and ghost propagators built into these calculations can 
be compared with Monte Carlo lattice simulations and are in excellent agreement lllj. 
While lattice calculations can only be performed with sizeable quark masses, the continuum 
Schwinger-Dyson/Bethe-Salpeter system can be computed even in the massless limit with all 
the essential physics of chiral logs built in. Consequently, this system provides a modelling 
of the chiral extrapolation ^3] so necessary to obtain physically meaningful results for light 
hadrons on the lattice. The successes described in Refs. 10-13 of this approach to strong 
physics justify the assumptions needed to truncate the Schwinger-Dyson equations and illus- 
trate how considerable progress has been made in extending the calculability of QCD from 
the perturbative regime to confinement scales so crucial for light hadron phenomena. 
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2 Scalar hadronic sector 



So far we have learnt that the dynamical breakdown of the chiral symmetry of QCD generated 
predominantly by uu , dd condensates ensures that pions are Goldstone bosons and the 
corresponding scalar field plays the role of the Higgs of the strongly interacting sector, its 
mass and that of all light hadrons reflecting the constituent mass of u and d quarks. But 
what is this scalar field? Is it the /o(400 - 1200) (or a), or /o(980), or /o(1370), or /o(1510), 
or /o(1720), or some mixture of all of these? None of these states is hkely to be a pure 
qq state, none likely to be pure glue, none solely qqqq or a KK molecule. All are mixtures 
of these, but what mixtures? This is the outstanding issue on which we try to shed a little 
light. 
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Figure 1: "Unquenching" of quark model states to make real hadrons has little effect on the 
vector mesons, p and </>, beyond allowing them to decay. In contrast the observed properties 
of the two scalar mesons /o(980) and ao(980) are produced by "dressing". These states then 
enhance the coupling of uu, dd systems to ss with no OZI suppression in scalar channels. 



The hadron states we know best, those that live the longest, are believed to be simply 
bound states of quarks. States of the quark model are most easily identified with the hadrons 
we observe experimentally when unquenching is unimportant. Fig. 1. Thus the (j) is readily 
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seen to be an ss state and the p and uo combinations of uu and dd. This follows from 
their respective decays to KK, and to 2ti and Svr. Though these decay modes are a crucial 
characteristic of their make-up, they have a relatively small effect on the states themselves. 
Thus only a small part of the Fock space decomposition of the physical is KK : it is 
predominantly ss. This is in part because of the P— wave nature of its hadronic dressing. 
The resulting small effect reproduces the suppression of the l/N^ expansion. In contrast, 
qq scalar mesons are strongly disturbed by their couplings to open hadron channels [T^ . 
Fig. 1. Thus almost regardless of their composition in the quenched approximation, the 
/o(980) and ao(980) are intimately tied to the opening of the KK threshold. Scalars change 
on unquenching. For them the l/N^ suppression of quark loops does not occur. The fact 
that these resonances, /o(980) and ao(980), couple to both vrvr/Tr?] and KK, means scalar 
non-strange and ss states communicate. Fig. 1. The coupling of different flavour quark pairs 
is not merely unsuppressed, nullifying the OZI rule in the scalar sector jT2j, but is even 
enhanced. This places strange quark pairs in the vacuum under the spotlight for further 
study dH. 

Since scalars are so intimately tied to the structure of the QCD vacuum, their nature is 
something we need to understand. Here we will continue to focus on the /o(980) and ao(980). 
It has been proposed for decades that these states have one of three possible compositions: 
either a simple qq structure (which for the /o(980) is dominated by an ss component, since 
we know it couples strongly to KK ), or a tightly bound four quark system or a looser 
KK molecule. A way to distinguish between these options is to study the /o and Qq in 
0-radiative decays. It has been advertised by Achasov [12], by Close and Isgur [20], and 
by others (see citations in Refs. 19,20) that these give rise to quite distinctive branching 
ratios, which for 7/o(980) are given in Table 1. There are analogous predictions for the 
ao(980): for instance, in the KK molecule picture, where K^K~" loops are key, the ratio 
BR{(j) 'yaQ)/BR{(j) 7/0) would clearly be one if the Qq and /o were degenerate in mass. 
We will see that the models used to predict the branching ratios for the 3 options shown in 
Table 1 are probably too simplistic, but that's for later. 

Table 1: Predictions for the absolute rate for —>■ 
7/0 (980) depending on the composition of the /o(980). 



Composition 


BR(0 ^ 7/o(980)) 


qqqq 


0(10-^) 


ss 


o(io-s) 


KK 


< 0(10-5) 



5 



(j) "y-KTv and — ^ 77177 have recently been measured at e~^e -colliders: at the Novosibirsk 
VEPP-2M with both the SND ^ and CMD-2 (2^ detectors and at Frascati with DAPHNE 
in the KLOE experiment The tttt and Tfrj spectra, Fig. 2, show a peaking at the end 

of phase space that the experiments identify with the /o(980) and ao(980), respectively. 
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Figure 2: 7r°7r° and tt^i] mass distributions in ^-radiative decays. The tttt results are from the 
KLOE PI, SND and CMD-2 detectors. The nr] data displayed are from KLOE 
in both the 77 and Sn decay modes of the rj. 



Let us look at how KLOE [zlj use the data of Fig. 2 to determine the decay rate of the 
/o(980) to compare with the predictions of Table 1. First they fit a Breit-Wigner form to 
the shape of the distribution, by adjusting the mass and width of the /q. Then one sees that 
this fails to describe the distribution at low vrvr masses, so part of the decay system must be 
provided by something other than the /o(980). There is, of course, the background to the 
same 'ymr final state from the sequential decays of — > pvr and p — > 77?. With the higher 
statistics from DAPHNE, the KLOE collaboration can separate out this particular decay 
by its distinctive angular distribution. One thus knows for the KLOE results that the 'ynir 
distribution (Fig. 2) is controlled by S-wave tttt interactions. As we have seen these cannot 
just come from /o(980) production. The remainder is taken to be the effect of the a. The 
cr's parameters are then taken from the Fermilab E791 experiment j2Sl with a mass of 478 
MeV and a width of 324 MeV. The amplitudes for a and /o(980) production are then simply 
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added and allowed to interfere. The resulting fit gives BR(0 7/o(980)) = (4.4 ±0.4) • 10""^ 
and one concludes the /o(980) is a four quark system — see Table 1. A similar analysis 
of the rjTT channel gives a branching ratio for the ao(980) in 0-radiative decay a factor of 6 
smaller. So why should we redo this analysis? 

Let us first concentrate on the isoscalar channel. The parameters of the /o(980) are not 
free variables. The same resonance pole position must appear in all processes to which the 
state couples. Consequently, one cannot permit the mass to be reduced by 10 MeV or more, 
or the width of a state that is typically found to be 40-60 MeV wide cannot be allowed to 
be 200-250 MeV wide. These changes made by KLOE increase the branching ratio through 
the /o(980) by an order of magnitude. To see why, consider the key dynamics of the process. 
The phase space involves a product of the vrvr momentum and the photon momentum. But 
because the photon is not any massless particle, but couples through a conserved current, the 
decay distribution in fact involves the cube of the photon momentum, as Achasov [^Hl has 
emphasised on many occasions. In terms of invariants this is proportional to (m|— M (vrvr)^)^. 
With the /o(980) so close to the end of phase space, small changes in its mass dramatically 
alter its branching ratio. This experiment cannot determine the parameters of resonances 
on its own. These parameters must however be the same as those required to describe 
other data. Moreover, the contribution of the a and the /o(980) must be added in a way 
consistent with unitarity. How to do this has recently been worked out by Elena Boglione 
and myself [77j . 

To proceed the basic assumption is that there are no strong (pre interactions. This means 
that the dominant strong interaction is between the final state pions. This is the presumption 
implicit in any isobar modelling of a decay. With the upper mass range fixed by the mass, 
there are very limited ways a tttt final state can be produced. Fig. 3. Either the radiates 
a photon leaving a tttt system that then interacts, or the (p radiates a photon producing a 
KK system that then interacts to produce a dipion pair. This occurs with the tt^tt^ system 
being in an isoscalar state, which is predominantly S'-wave. Then coupled channel unitarity 
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Figure 3: The unitarity constraint relevant to tttt interactions in a definite partial wave in 
0-radiative decay. When there are no significant (pn strong interactions, or these have been 
removed from the data, the constraint is exact. Its solution is Eq. (1). 
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Figure 4: Plots (a) and (b) show the I = J = ReVAMP hadronic amphtudes j2H] 
|T(7r7r 7T7i)\ (labelled Tn) and \T{KK tttt)! (labelled T21) from which F{s) is con- 
structed, according to Eq. (1), where ^/s = M^tttt), the cm. energy. 

requires that the amphtude, F for ^ •jtctt can be related to the basic amplitudes T for 
TTTT TTTT (Tn) and KK — s> TTTT (T21) by 

F{(f) —>■ 77r7r) = ai{s) T(7r7r tttt) + a2{s) T{KK —>■ tttt) , (l) 

where s = M(7r7r)^. With no strong 07r interactions, the two functions aj(s) are real. They 
can be interpreted as the intrinsic couplings for ^-k-k and — '-/KK, respectively for 
z = 1,2 (Fig. 3). Since the whole point of studying this process is because the initial state 
is almost 100% ss , we expect the coupling 0^2 to be much larger than ai. Our analysis 
shows that experiment does indeed support this. However, the functions ai{s) each have 
a factor of (m| — s), as previously explained, following from QED gauge invariance. This 
strongly reduces the contribution from T{KK tttt), which is dominated by the /o(980), 
and enhances the contribution from T(7r7r tttt), which is controlled by the /o(400 — 1200) 
(or a), and from which the /o(980) effectively decouples — see Fig. 4. So while there is a 
sizeable /o(980) component, much of the nn decay distribution is produced by tttt interactions 
outside the narrow /o(980) region. 

By building into the analysis known experimental information on the scattering reactions 
TTTT TTTT and TTTT —>■ KK, which automatically embodies details of the /o(400 — 1200) and 
/o(980), we can use the data on — > 77r°7r° to determine the couplings of these scalar 
resonances to this channel in as model-independent a way as possible. This is the purpose 
of the recent analysis by Elena Boglione and myself [27J. As a simple template we first 
used the old hadronic amplitudes determined by David Morgan and I j,28j (called ReVAMP 
as explained in Ref. 27). These have the /o-pole on sheet II at M(7r7r) = (988 — i ■ 23) 
MeV. Factoring out the Adler zero and the photon momentum required by QED gauge 
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invariance for the radiative decay process, we then have constant couphng functions, ai{s), 
and obtain the fit shown in Fig. 5. The quahty of the fit is excellent indicating no reason 
to expect significant strong (pre interactions need be included and showing the final state vrvr 
interactions in this decay are completely consistent with those from other processes built 
into the Re VAMP amplitudes. 




I I I I I I I I I I 
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E 

Figure 5: Simultaneous fits to the data on the tt^tt" decay distribution in — > 77r''7r'^ 
from the KLOE |^ and SND |^ collaborations. These fits have been obtained using the 
Re VAMP set of underlying amplitudes, with just 3 real parameters ^7\, as described in the 
text. E = M{ti-k). 

It should be clear that the exact branching fraction for the /o(980) is exceedingly sensitive 
to the position of the corresponding pole, because of its nearness to the edge of phase space. 
Moreover, one does not see the /o(980) as a simple Breit-Wigner shape (compare Fig. 4b 
with Fig. 5), so in fact its branching fraction is not directly related to an experimental 
observable. The only well-defined quantity is the residue at the pole which truly represents 
the coupling — >■ 7/9. For a state that overlaps with another broad resonance and with a 
strongly coupled threshold, the idea of a branching fraction is wholly model-dependent. 



9 



The fit shown in Fig. 5 gives a — 7/0 (980) coupling of 0.62 MeV (see Ref. 27 for 
the exact definition), which approximates to a branching ratio of 10% of the total tctt S- 
wave decay being through the /o(980), i.e. a branching fraction of 0.31 ■ 10"'^. Looking at 
Table 1 we see this corresponds no longer to the qqqq composition inferred by KLOE 
but being an order of magnitude smaller is like that for an ss constitution. However, this is 
too simplistic. Oiler ^] has noted that the KK molecular picture can include not just the 
charged kaon loops used in the predictions of Table 1, but also neutral kaons too. Then the 
prediction extends from 10~^ to 10~^ to encompass our result too. 

What of (f) 7ao(980) '-/{tt^t])'! A coupled channel analysis like that for the tt^tt^ 
decay mode is not presently practicable because of lack of precision information on ni] ttt] 
and TTT] KK scattering in S'-wave channels. However, these have none of the complication 
of overlapping resonances of the isoscalar tctt mode. Consequently, it is more likely that 
such an analysis when possible will reveal a BR{(f) 709 (980)) much more similar to that 
found by KLOE j2Sl- This is the same order of magnitude as our smaller BR{(f) 7/0 (980), 
which is something expected if there is a large KK component of these two scalar mesons. Of 
course, the most simple minded interpretation of the KK molecule picture with just charged 
kaon loops seeding the decay gives the ratio of branching ratios for oq and /o to be one, only 
if these states are degenerate in mass. Just a 10 MeV difference in their masses changes this 
to 0.4 or 2, depending on the sign of the mass difference. Consequently, model predictions 
and analyses of data are extremely sensitive to the /o(980) and ao(980) pole positions. 

Elena Boglione and I have illustrated this by considering in our analysis of the 7r°7r° 
data the more recently determined hadronic amplitudes of Anisovich and Sarantsev ^Uj . 
These give the /o-pole at M(7r7r) = (1024 - i ■ A3) MeV. Fits displayed in Ref. 27 show that 
this much greater width for the /o(980) results in a coupling of 1.9 MeV, which translates 
into a modelled branching ratio of 1.9 ■ 10~^ — a factor 6 larger than using the ReVAMP 
amplitudes with a much narrower /q described earlier. However, the use of Anisovich and 
Sarantsev amplitudes (AS) gives fits of much poorer quality. If it were not for the fact that 
these authors treat a far greater range of more recent data, such as that from Crystal Barrel 
and from GAMS (see Ref. 30 for details), than in the older ReVAMP analysis ^Ei, the 
quality of fit would dismiss such a wide /o(980) as unlikely. Nevertheless, the fact that the 
7/0 coupling is so sensitive to the details of the /o-pole position means that we need 
to tighten up our determination of the underlying hadronic amplitudes in this crucial mass 
range to be certain of the results. This would not be the case if this same final state were 
to be studied in the decay of the recurrence of the (p at 1680 MeV. This may indeed become 
possible in photoproduction in Hall D at Jefferson Lab. 

States in the spectrum are identified as poles of the ^-matrix. Crucially, their position in 
the complex energy plane is independent of the process in which they appear. Experiment 
does not, of course, deliver direct information about such complex poles, but only about 
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real quantities on the real energy axis. It has been argued [3T] that this means poles are 
irrelevant. Only experimental observables matter and these will eventually be calculable from 
QCD on the lattice. Such calculations may indeed reproduce experiment precisely. However, 
agreement does not necessarily bring understanding. Monte Carlo lattice simulations are a 
black box. A perfect prediction for the e'^e~ tt+tt^ cross-section and the J = 1 vrvr 
phase-shift below 1 GeV gives no comprehension that the p-resonance with a universal pole 
position describes both datasets. More starkly, only the universality of poles tells us that 
the appearance of the dip near 1 GeV in J = J = tttt scattering (Fig. 4a) and the peak 
in ^ 7(7r7r) (Fig. 2) are manifestations of the same /o(980). However accurately lattice 
QCD computations reproduce data, only continuation into the complex energy plane brings 
understanding of the spectrum of hadrons and their couplings. Consequently, precision 
knowledge of the the pole positions and residues of the /o/ao(980) is essential if we are to draw 
conclusions about the composition of these states in as model-independent way as possible 
— independently of the modelling of 0-radiative decay by Achasov JH] or the generation of 
these states by unitarisations of chiral dynamics [3211211 — just using experiment. 

A rich source of additional information on light hadron final states is being provided 
by the decays of heavy flavour mesons. That final state interactions and the way they 
universally appear shapes such decays has recently been studied in D — Svr [33J- Data are 
typically analysed in an isobar picture |25l , where it is assumed the three final particles 
interact only in pairs with the third as a spectator. If three body forces are needed, they are 
usually assumed to have constant matrix elements and so populate the Dalitz plot according 
to phase-space. In the analysis of the Fermilab E791 results |2S1 with ~ 1100 —>■ tt+tt+tt^ 
events, the known resonances that couple to tttt — the p, /2(1270), etc. — are included. 
The resulting fit is poor at low tttt masses. However, this is dramatically improved if an 
I = J = resonance of mass (478 ± 24 ± 17) MeV and width (324 ± 42 ± 21) MeV is added 
in Breit-Wigner form. Hence the E791 group claim to have confirmed the a resonance 
This is to forget that, assuming an isobar model, one has by definition now determined the 
/ = J = TTTT — s> TTTT iutcraction. Though the phase-shift is large, it is undoubtedly not given 
alone by a simple Breit-Wigner form with the claimed mass and width. 

With ~ 1500 events on this same decay, the FOCUS group 33 at Fermilab has con- 
firmed the E791 data and analysis. But if instead the low mass S'-wave vrvr interaction is 
parametrised using the Anisovich and Sarantsev description, then the fit to the D 3tt 
Dalitz plot is even better. Though here and in Ref. 27, we have queried the wide /o(980) 
of AS, this matters little for D —>■ Sn decay. Their amplitudes and those of Re VAMP are 
virtually the same from vrvr threshold to 900 MeV (see Fig. 1 of Ref. 27). What Malvezzi 
and collaborators [321 have shown is that the vrvr final state interactions in D ^ 3vr decay 
are completely consistent with what we know of such interactions from all other processes. 
They embody nothing new! 
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Whether there is a cr or not as a short-hved resonance is then not a question of whether 
it is just seen in D ^ 3tt decay, but whether it is there in all other I = J = tttt final states 
too. The difficulty in answering the question of whether data on the real axis have sufficient 
precision to determine the existence of such a very distant pole jHlI is amply illustrated by 
the difference between the AS [30] and Re VAMP j2Hl amphtudes. AS have no low mass pole, 
while ReVAMP does : yet each describes essentially the same experimental results below 
900 MeV. 

What we learn is that analysing data on a single channel, where ffnal state interactions 
are important, cannot be meaningfully done in isolation. Unitarity requires consistency 
between reactions. Only by analysing data from different processes with the same final states 
simultaneously can we hope to be able to draw definitive conclusions about the fascinating 
scalar sector. Since these states with zero quantum numbers refiect the nature of the QCD 
vacuum, further study is essential. Precision results from Hall D at Jefferson Lab should 
contribute enormously to this endeavour. 
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